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Abstract 

We study effects of fluctuations on the mesoscopic length-scale on systems with mesoscopic inho¬ 
mogeneities. Equations for the correlation fnnction and for the average volnme fraction are derived 
in the self-consistent Gaussian approximation. The equations are further simplified by postulating 
the expression for the strnctnre factor consistent with scattering experiments for self-assembling 
systems. Predictions of the approximate theory are verified by a comparison with the exact resnlts 
obtained earlier for the one-dimensional lattice model with first-neighbour attraction and third- 
neighbour repulsion. We find qualitative agreement for the correlation function, the equation of 
state and the dependence of the chemical potential y on the volume fraction C Our results conhrm 
also that strong inhomogeneities in the disordered phase are found only in the case of strong re- 
pnlsion. The inhomogeneities are reflected in an oscillatory decay of the correlation function with 
a very large correlation length, three inflection points in the /r(C) curve and a compressibility that 
for increasing ( takes very large, very small and again very large valnes. 
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I. INTRODUCTION 


Density functional theory (DFT)[l| is very successful in describing microscopic properties 
of fluids. However, when fluctuations at length scales much larger than the molecular size 
dominate, predictions of the standard DFT are much less accurate. Fluctuations at meso¬ 


scopic lejigth scales are present in a critical region , at an interface between coexisting fluid 


phases [3| and in systems with conmeting interactions. Competing interactions are quite 
common in soft matter systems jd-^. For example, charged particles repel each other at 
large separations, but at short separations attract each other with solvent-mediated forces. 
The above so called SALR potential may lead to self-assembly into clusters of spherical or 


elongated s^ 


m space 


molecules 


rape, into layers or into networks, and to a periodic distribution of these objects 


12 


13 


. Competing interactions are also present in systems containing amphiphilic 
14| . and in thin magnetic hlms with dipolar interactions 1^, Il6 |. 

For systems with competing interactions mean-held (MF) theories, including a local 


DFT, predict stability of severa 


temperature phase diagram 


0 


ordered phases inside a bell-shaped region in the density- 


17 


20|. The disordered phase stable outside this region 


coexists with two ordered phases only - with periodically distributed clusters at low- and 
with periodically distributed voids at high density. The coexistence line between the dis¬ 
ordered and the ordered phases lies close to the line of instability of the disordered phase 
with respect to periodic density huctuations. When the instability line is approached, both 
the correlation length and the amplitude of the correlation function diverge HQ. The 
latter result is clearly nonphysical, since the correlation function for dimensionless density 
is bounded from above by 1. 


In simulations the ore 
signihcantly different 


ered 


22 


phases have been observed too, but the phase diagram is 


23| . The stability region of the ordered phases is much 


smaller [l2, 22, 23|. In addition, a less ordered phase is stable between the stability regions 
of two different ordered phases Q,H|, in contrast to the MF results Q,[u, Q- In the case 
of a triangular lattice model with competing hrst-neighbor attraction and third-neighbor re¬ 
pulsion at low temperature T a molten lamella phase occupies some phase-space region 
between the stability regions of the phase with periodically distributed clusters and the phase 


with periodically distributed stripes 


22| . At higher T the disordered fluid is stable between 


the ordered cluster phase and the molten lamella. The transition between the fluid and the 
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molten lamella is continuous at high T and becomes first order at the tricritical point. In 


the molten lamella discovered in Ref. 


22| the orientational order of stripes is present but 


the translational order is absent. A similar phase called a “nematic phase” was recently 
detected in magnetic films with competing interactions 2^ in a continuous Brazovskii-type 


model 


25l |. According to these new discoveries, instead of the fluctuation-induced first-order 


phase transition between the disordered and the periodically ordered phases predicted ear¬ 


lier by Brazovskii 


25| . a continuous or a first-order transition between the fluid and the 


molten-lamella phases occurs. 

In addition to the presence of new phases and signihcant modihcations of the phase- 
coexistence, the fluctuations at the mesoscopic length scale lead to significantly different 
properties of the disordered inhomogeneous phase. In simulation snapshots of the disor¬ 
dered phase clusters or layers are clearly visible, indicating that this phase can be very 


inhomogeneous j8|, l9|, l22|, [2^. Even though the distribution of these objects is not peri¬ 


odic, they can be correlated over large distances. The disordered inhomogeneous phase may 
thus be considered as a molten periodic phase, with the long-range order destroyed by the 
mesoscopic fluctuations. 

The effects of fluctuations on the disordered phase were determined in a one-dimensional 
(Id) lattice model with competing first-neighbor attraction and third-neighbor repulsion by 


the exact transfer matrix method 


18[ |. In the ground state (T = 0) of this model a dilute 


phase (vacuum) is stable for small values of the chemical potential /i and a dense phase 
(fully occupied lattice) is stable for large values of p. When the repulsion is sufficiently 
strong, an ordered phase with periodic distribution of clusters and density p = 1/2 is stable 
for intermediate values of p. At T > 0 only a disordered phase is stable, but the presence 
of clusters at low T is reflected in a very large correlation length, and in a specific shape 
of the p{fi) curve. It contains three inflection points at low T - the central one at p = 1/2 
and the other two at p corresponding to the T = 0 phase transitions. For T —)■ 0 a step- 
like p(p) curve was obtained. The very small compressibility of the system for p 1/2 
indicates formation of clusters that repel each other at short separations. The very large 
compressibility at the other two inflection points signal an approach to the phase transitions 
that occur at T = 0. Only at high T a single inflection point at p = 1/2 is present. In 
contrast, in ME the p(p) curve has a single inflection point for the whole range of stability 
of the disordered phase. Thus, ME fails to predict the qualitative features of the disordered 
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inhomogeneous phase. Only for high T, where strong inhomogeneities are no longer present, 
qualitative properties of the disordered phase are correctly described by MF. MF predicts 
stability of a periodically ordered phase instead of strongly inhomogeneous disordered phase. 

From the comparison of the MF and simulation or exact results it evidently follows that 
an inclusion of the most relevant mesoscopic fluctuations is necessary not only for the quan¬ 
titative, but also for the qualitative description of systems with mesoscopic inhomogeneities. 
A tractable theory that would accurately predict properties of systems with inhomogeneities 
on the mesoscopic length scale is still missing, however. A theory correctly incorporating 
mesoscopic fluctuations should predict at least a correct topology of the phase diagram and 
correct properties of the disordered phase. 

Effects of periodic fluctuations of the order parameter (OP) 0 are described by the Bra- 
zovskii held theory j^. However, the Brazovskii functional F_b[ 0] depends on free parame¬ 
ters, and quantities such as compressibility for a particular volume fraction ( and T cannot 
be determined, unless a relation between the free parameters and measurable quantities is 
known. An attempt to express the free parameters in L r [d)1 in terms of T and ( was made 
in Ref.[l4, 17, 26|. The theory developed in Ref. jl4 . 17 . 261 is rather complex, however, and 
additional approximations are necessary to solve the obtained equations. 

The purpose of this work is further development of the mesoscopic description of in- 
homogeneous systems that combines the DFT and held-theoretic approaches 


13 , 


27|. 


We introduce a tractable approximation for the disordered phase. In sec.2 we present a 
derivation of a self-consistent equation for the correlation function, and the equation for 
the average volume fraction in terms of T and /i. The equations are obtained by the DFT 
methods, and the contributions associated with the mesoscopic fluctuations can be calcu¬ 
lated by the methods of held-theory. We show that the term oc 0^, usually neglected in 
the Brazovskii-type theories, plays a very important role in fluids and soft-matter systems. 
In sec. 3 we derive the equations in the self-consistent Gaussian approximation. In sec. 
4 we limit ourselves to the disordered phase, make further assumptions and obtain much 
simpler equations that can be solved easily. In sec. 5 our equations are solved for a Id model 
with the SALR potential. The results of our theory are compared with the exact solutions 
obtained in Ref. H for a Id lattice model. We discuss the accuracy of our approximations 
in sec. 6. 
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II. DERIVATION OF THE MESOSCOPIC DENSITY FUNCTIONAL THEORY 


We consider systems with competing interactions, where mesoscopic inhomogeneities oc 
cur on a length scale A significantly larger than the size of molecules a. Following ref. 17l.l26| 
we collect all the microscopic states into disjoint sets; each set represents one mesostate de¬ 
scribed by a smooth function C(r)- C(i') is equal to the fraction of the volume of a sphere 
with a center at r and a radius A S> i? ^ a that is covered by the particles in each microstate 
belonging to the set represented by C(i')- Another words, the mesoscopic state represents all 
the microscopic states that are obtained by changing the positions of the particles without 
changing the volume occupied by the particles in each mesoscopic region. By fixing the 
mesostate (’(r) we impose a constraint on the available microstates. The grand thermody¬ 
namic potential in the presence of this constraint is denoted by Uco[C(r)]- exp(—/Sr2co[C(i')]) 
is equal to the sum of the Boltzmann factor over all microscopic states with frozen 

mesoscopic fluctuations. In the above 13 = l/lksT), ks, T and H are the Boltzmann fac¬ 
tor, temperature and the Hamiltonian respectively. The grand potential in the presence of 
mesoscopic fluctuations is given by 


U = —/coTln 


( 1 ) 


where 

S = y (2) 

and the functional integral is over all mesoscopic states that by definition are C < 1- Since 
the summation of over all microstates compatible with (j'(r) is included in and 

in ([2]) we perform the summation over all mesostates C(i ')5 ^ contains contributions from 
both, microscale and mesoscale fluctuations. 

We introduce mesoscopic fluctuation by 


0:=C-C 


(3) 


where C, denotes the average volume fraction, and rewrite ([T]) in the equivalent form 


with 


H = Ueo[C] - ksTln ^ j 

(4) 

Hj[C4>] = n^[C + 4>]-n^[C]. 

(5) 
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Eq. (jl]) contains two contributions, and each of them is associated with the fluctuations on 
a different length scale. In the hrst term the mesoscopic volume fraction has its equilibrium 
form. This term contains contribution from the fluctuations on the microscopic length scale. 
The second term contains the contributions from the fluctuations on the mesoscopic length 
scale, i.e. from different mesoscopic inhomogeneities that are thermally excited with the 
probability /S. When ( is the average volume fraction, then it follows that (0)t = 0, 
where 

/ Dc^Xe-f^^f 

f 

In the following = 0 is always assumed. 

Our purpose here is the development of an approximate theory that allows for a determi¬ 
nation of ( for given temperature and chemical potential with the fluctuation contribution 
in (jlD taken into account. We shall also hnd approximations for the equation of state (EOS), 
the pair correlation function and the boundary of stability of the disordered phase. Note 
that both terms on the RHS of Eq. (Hj) are functionals of (. The equilibrium volume fraction 
( corresponds to the minimum of the grand potential i.e. in equilibrium the hrst functional 
derivative of 0[C] w.r.t. ( must vanish, and the second functional derivative must be positive 
dehnite. We introduce the functional derivatives 





Cn{ri, ...,r„) := 

In 

5C(ri)-(5C(rn) 


_ 

SC{ri)...SC{rn) 


where dH) was used, 

<5C(ri)...5C(rn)’ 

and we do not indicate the functional dependence of and Cn^ on (. 
For n = 1, 2 we obtain from ([7|) 



and 


where 


C'li 




C2(ri,r2) = C^°^(ri,r2) ( 




'Jf ~ XWTTTWTTT// 


5C(ri)(5C(r2) 5C(ri) (5C(r2)' 


{XYYp := {XY)f - {X)f{Y)f. 


(7) 


( 8 ) 


(9) 


( 10 ) 

( 11 ) 


6 









We neglect the fluctuation contribution to Cn for n > 3, i.e. we make the approximation 

Cn = C® for n > 3. (12) 

From (|9]) and the requirement Ci(r) = 0 we can obtain ^ for given T and the chemical 
potential fi, if we know the form of flco and we can perform the functional integrals in the 
second term on the RHS of ([9]). 

We assume the standard focal mean-held approximation for the grand potential with 
suppressed mesoscopic fluctuations, 


«o.[Cl = C'IC]-rS|c]-fiiV|c]. 


(13) 


The entropy S in the local density approximation for fixed mesoscopic volume fraction is 
given by 

-rS[cl = |*A(C(r)). ( 14 ) 

Different approximations for the free-energy density of the reference (hard sphere) system, 
fh{C), can be chosen. The internal energy for fixed mesoscopic volume fraction is given by 

fo[C] = ^ 

where V (ri — r 2 ) = M(ri — r 2 )g{ri — r 2 ), with u and g denoting the interaction potential and 
the pair distribution function respectively. Note that we use volume fraction rather than 
density in flTHll . therefore we should re-scale the interaction potential M(ri — r 2 ) by the factor 
(b/vr)^ to obtain the same energy as in the standard theory. We also re-scale the chemical 
potential, /2 = (6/7r)/i, so that 

h^[C] =f^J drCir), (16) 

where N[(] is the number of particles for given (. In this local density approximation we 
have 

CfV) = J driC(ri)/3I/(ri - r) -F Ri(C(r)) - l3fi (17) 

C'f^(ri,r2) = /3fo(ri - r2) R2(C(ri))(5(ri - r2) (18) 

and for n > 3 

•••,rn) = A„(C(ri))(5(ri - r2)...(5(r„_i - r„), (19) 


where 


7ln(C) 


dC^ 


( 20 ) 
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In order to perform the functional integrals in (|1]), dH]) and (fTU|) we expand Hf[(,(f)] in a 
functional Taylor series w.r.t. 0, 

fiHf[C,cl>]=fiHo + fiAH ( 21 ) 

with 

l3Ho[C,(p] = ^ j dvi j dr20(ri)C'f Vn r2)0(r2) (22) 

and 

/5Aff[C,0]= f drci°\r)(p{r) + Y^ j (23) 

We truncate the expansion in 0231) at the fourth order term. This truncation is justified, 
because by definition the volume fraction is less than 1 , and for small fluctuations, 0 <C 1 , 
the higher order terms in fl2^ are irrelevant. The functional integral in (jl]) can be extended 
to arbitrarily large functions 0 , because the large fluctuations are strongly damped in (j4]) 
by the Boltzmann factor for Hf given in fl2Tl) - fl2^ , and do not influence the results 

in a signihcant way. 

Note that when Hj is given by fl?TD - fl2^ , the fluctuation contributions in ([2]) and flTU]) 
can be expressed in terms of (0"(r)) and ( 0 ”(ri) 0 ™(r 2 )). Thus, it is necessary to calculate 
the correlation functions for the mesoscopic fluctuations of the volume fraction in order to 
obtain the expressions for Ci and C 2 . In order to calculate these functions field-theoretic 
methods could be used. In the next section we present an approximation that allows to 
determine C 2 and ((/)(ri) 0 (r 2 )) in a relatively simple way. 

III. SELF-CONSISTENT GAUSSIAN APPROXIMATION 

In practice we can calculate the Gaussian functional integrals only. We approximate Hf 
by an effective functional that is quadratic in the fluctuation 0 , 

^Hf ^ (3Hg = ^ j dri j (ir20(ri)G(ri,r2)0(r2). (24) 

In this approximation the form of C differs from C^\ so that the effect of AH (see 
is indirectly included in the average quantities (see ([ 6 ])). Note that when Hf is approximated 
by Hg, then the correlation function G(ri,r 2 ) := {(j){ri)(j){r 2 )) is given by 


dr2G(ri,r2)G(r2,r3) = (5(ri,r3). 


( 25 ) 



The best approximation for C(ri,r 2 ) is such that G in is as close as possible to the 


exact correlation function, 
inverse correlation function 


^his means that C should be as close as possible to the exact 
ij. Here, we assume that C = C 2 - Note that in the Gaussian 
approximation the 2 n-point correlation functions can be expressed in terms of products of 
the two-point correlation functions. Thus, because the fluctuation contribution in Eq. ffTOj) 
consists of terms proportional to ( 0 "'(ri) 0 ™'(r 2 )), it relates C 2 = C with G, and Eq. fl2H]) 
relates G with G. In this self-consistent Gaussian approximation (X) j in Eq. fflU]) is replaced 
by {X) = J D(pX exp{—(3HG)/ f D(/)exp(—/3HG), with Hq given in (]2T|) . 

In order to calculate the functional integrals in Eqs.([9]) and (1T0|) in the above Gaussian 
approximation, we note that from fl 2 T|l . ([7]), ffT 2 |) and flT^ it follows that 


It = 


(26) 


and 




rj). (27) 

and oa, using flTTI) - fl 2 ni) and performing the 


After inserting the above equations in 1 
Gaussian integrals, we obtain the main results of the self-consistent Gaussian approximation 


fj/x = f driC(Ti)fiV(vi - r) + 2 l,(C(r)) + r) 


(28) 


and 

C(ri,r2) = Gt(ri,r2)-F 


(o)^„ „ ^ ^ "43(C(ri))A3(C(r2))^2 


-G(ri,r2)(5(ri -r2) 


G\ri,r2). (29) 


2 V V ^ 2 

Eqs. fl29D and fl25D have to be solved self-consistently. The last term in fl2^ and the 
last two terms in fl2^ represent the contributions from the mesoscopic fluctuations. We 
should stress that because of the coarse graining, the correlation function for the mesoscopic 
volume fraction, G(ri,r 2 ), differs from the correlation function for the microscopic density. 
In Ref.ln, 26| it was shown that G(ri,r 2 ) represents the microscopic correlations averaged 
over mesoscopic regions (smaller than the scale of inhomogeneities) around the points ri 
and r 2 . As a consequence, G(ri,r 2 ) 7 ^ 0 for ri = r 2 , because it contains contributions from 
the microscopic correlations between the particles with the centers belonging to the sphere 
of the radius R larger than a/2, i.e. between the particles which do not overlap. Note also 
that the fluctuation contribution leads to a shift of the average volume fraction for given T 
and ft compared to the ME prediction. 
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FIG. 1: The Feynman diagrams representing the fluctuation contribution in Eq. (|28|) (a) and in 
Ea. (l29p (b) and (c). Thick line connecting vertices at ri and r 2 represents G(ri,r 2 ). The vertex 
with a triangle or the square represents A 3 or A 4 respectively. 


Our result fl2^ agrees with the well known expression in the field theory [2^ at the 
self-consistent one-loop approximation. The corresponding Feynman diagrams are shown in 
Fig.l. In the standard field-theoretic approaches, however, are free model parameters, 
and the last term in fl2^ is usually neglected. In our theory A 3 is a function of C and 
vanishes only at a single value of the volume fraction, corresponding to the critical point in 
a system with attractive interactions. In our theory ^3 plays a very important role for the 
relation between the average volume fraction and the chemical potential. 


IV. SELF-CONSISTENT GAUSSIAN APPROXIMATION FOR THE DISOR¬ 
DERED PHASE 

In this section we restrict our considerations to a disordered isotropic phase with the 
mesoscopic volume fraction independent of the position, and to isotropic interactions, i.e. 
U(r 1 — r 2 ) depending only on r = | ri — r 2 1 . In this phase G and C depend only on r = | r 1 —r 2 1, 
and we simplify the notation, introducing G{r) = G(ri,r 2 ) and G{r) = C'(ri,r 2 ). Due to 
the translational invariance, Eq. 02811 takes the form 

DP = C j drDVir) + A.iO + ^5, (30) 
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and we can write Eg. 029 p in Fonrier representation 



(31) 


where F denotes the Fourier transform of F, and we have introduced 


g = G(r, r) 


(32) 


and 


D{k) 


j dvG^{ry^-\ 


(33) 


In the disordered isotropic phase Eq. fl25p takes the form 


C{k)G{k) = 1. 


(34) 


In the context of fluids with competing interactions, ^ 0 except from the critical 
value of the volume fraction, y ^ 3 (C) < 0 for (^ < (c, and A 3 (C) > 0 for ^ > (c- Thus, 


the fluctuation contribution in fl30p leads to decreased and increased value of the chemical 


potential for C < Cc and C > Cc respectively compared to the MF prediction, provided that 
Q > 0. This effect of mesoscopic fluctuations is independent of the details of the interaction 
potential. Note that the value of Q is the larger the stronger are the inhomogeneities on the 
mesoscopic length scale. 

In order to solve Eq. flTip for a particular system, we need to know the form of the 
interaction potential. However, we can make general qualitative or semi-quantitative pre¬ 
dictions for a class of systems with interactions such that V{k) has a pronounced minimum 
at k = ko > 0. The wavenumbers of the dominant mesoscopic fluctuations correspond to 
the maximum of the Boltzmann factor exp(—/SFc), and we assume that onlv k ^ kn are 


relevant when the peak in exp{—f3V{k)) is pronounced and narrow 
such systems, and make the approximation 



Ci°\k) ss Co + Vaik^ - k^f 


(35) 


where 


Co — A 2 + PV (fco) 


(36) 


and 


/9F"(A)o) 


(37) 
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We have used (fTSD and have made the approximation fc + /cq ~ 2/co, valid in the neighbor¬ 
hood of ko- The above assumption is justified for an even function of k in the neighborhood 
of ko > 0, where the first derivative w.r.t. k is = 0, and the second derivative is 

l/"(fco) > 0. The results obtained with this assumption concern any system with inhomo¬ 
geneities on the mesoscopic length scale ~ 271/ko strongly favoured energetically compared 
to inhomogeneities on different length scales. Thus, the results presented below concern 
systems with a deep minimum of V{k). For potentials with a shallow minimum the approx¬ 
imation (155]) can be an oversimplihcation. 

The equations fl3Tl) - fl34D are still rather difficult. It is thus reasonable to make further 
approximations in order to verify if the predictions of this theory are correct at least on a 
qualitative level. Instead of solving numerically Eqs. (l3TD - (l34|) . we postulate a particular 
form of C{k), and determine the parameters in the postulated expression. Experimental 
results for the structure factor in systems with mesoscopic inhomogeneities are reasonably 
well described by a function inversely proportional to C^\k) given in Eq. (I35|) . with co,ko 
and vq treated as fitting parameters. Based on this observation we assume for C(k) the 
same form as in Eq. (l35|l . but with renormalized parameters, i.e. we postulate 

C(k) TH Cr + Vr(k^ — k^y. (38) 


In order to determine the renormalized parameters we need 3 equations for the 3 un¬ 
knowns. We require that the first derivative of the function given by Eq. fl3Tl) vanishes at 
k = kr- The value of this function at its minimum is c^. Finally, we require that the second 
derivatives at kr of both expressions, Eqs. fl31l) and fl55]) . are the same. The three require¬ 
ments guarantee that the shape of C(/c), Eq. fldTl) . is reproduced by Eq. fl38l) at least near the 
minimum. From these requirements we obtain the set of equations 


Cr = Co+ Vo{kr - k^f + - ^D{kr) 

Avoikl - kl)kr = ^D'{kr) 

Svrkl = 4vo{3kl - kl) - ^D'\kr). 


(39) 

(40) 

(41) 


Note that in the standard Brazovskii approximation, i.e. with = 0, from (1401) and 
(ITT]) we obtain kr = ko, Vr = vq, and Eq. (l39]) can be easily solved analytically j^. 
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The pressure p = —fi/V, where V is the system volume, in this theory is given by (see 

(HD, m-m) 

P = -iv'(0)C" - A(0 +^(+^''^[1 D4>e-^’'A . (42) 


We express p in terms of ( according to Eq. 
from (I12D the EOS of the form 


, evaluate the functional integral and obtain 


P = - A(C) + ksTA.iOi + ksT 


'A.iocg 

/ In 

(CikR 

2 

1 2{2tiY 

1 27r / 


(43) 


The integral in fl43D is over the spectrum of mesoscopic fluctuations, and is cutoff- 
dependent for C{k) given in Eq. flSSp . In the mesoscopic theory the cutoff is naturally pro¬ 
vided by the scale R of coarse-graining. However, the pressure should be independent of the 
coarse-graining procedure and the cutoff-dependent part should cancel against the neglected 
cutoff-dependent contribution to hlco- We find that the cutoff-independent contribution to 
the integral in fl43D is proportional to ^Cr/vr for the 3d and Id systems. Note that the 
approximations fl35D and fl38D are valid for strong inhomogeneities at a well-defined length 
scale, i.e. for Cr/vr 1. Since Q oc [2^ (see also the next section), in the range of 

validity of the approximations, the last term in fHSp is negligible compared to the fluctuation 
contribution proportional to and will be disregarded. 


V. RESULTS FOR THE SALR POTENTIAL IN ID 


In this section we test the accuracy of the self-consistent Gaussian approximation with 
the further approximations described in sec. 4 by comparing our results with the exact 


solutions obtained in Ref. 


18| for a Id lattice model with competing attractive and repulsive 


interactions between the first and the third neighbors respectively. In the lattice model 
solved exactly in Id the ratio between the third-neighbor repulsion J 2 and the hrst-neighbor 
attraction Ji is denoted by J, and the interaction potential in Fourier representation takes 
the form 


f3V{k) = 2(3* ( J cos(3/c) — cos(/c) 


(44) 


where (3* = 1/T*, and T* = ksT/ Ji is temperature in units of the hrst-neighbor attraction. 
V{k) in Eq. 0441) assumes a negative minimum for kQ > 0 if J > 1/9. By V*{k), p* and p* 
we shall denote the corresponding quantity in units of R. 
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From the form of fEqs. fl5^ - fl57|) ) it immediately follows that in MF the disordered 
phase is unstable for T* < —V*{ko)/A 2 {C). The explicit form of ^2(0 follows from (1201) and 
the reference-system free energy of the lattice-gas form (note that in Id the volume fraction 
and the dimensionless number density are identical), 


/3A(C) = ClnC + (l-C)ln(l-0- 


(45) 


This well known MF result 


HHQ is i 


incorrect, since in Id models with short-range 


interactions there are no phase transitions for T > 0 


28j. Exact results |l8| indicate. 


however that the disordered phase is strongly inhomogeneous in the phase-space region 
where MF predicts its instability. Thus, it is interesting to verify predictions of our theory 
for T* < -V*{ko)/A 2 {C). We shall focus on T* < -V*{ko)/A2{0 for J = 3 and J = 1/4, 
because most of the exact results of Ref. concern J = 3 and J = 1/4. The phase-space 
region of interest is enclosed by the solid and dashed lines in Fig.2 respectively. 

The two chosen values of J correspond to systems with qualitatively different properties. 
In the case of J = 1/4 only the empty or the fully occupied lattice is present in the ground 
state (GS) for small or large values of p* respectively. In a system with the short-range 
attraction stronger than the long-range repulsion, exact results show stability of the 
disordered phase with oscillatory and monotonic decay of correlations at large and at very 
small T* respectively. The compressibility at ~ 1/2 is very large at low T*, signaling the 
approach to the phase separation at T* = 0. 

For strong repulsion (J > 1) clusters composed of 3 particles separated by 3 empty sites 
are favoured energetically, since there are as many attractions (occupied nearest-neighbour 
sites) as possible in the absence of repulsion. Such a periodic phase is stable for intermedi¬ 
ate values of fi* in the GS. Exact results show that for J = 3 the correlation function 
exhibits oscillatory decay with the wavelength ~ 6. The correlation length is several orders 
of magnitude larger than the particle size for ( cs 1/2 and T* —R(fco)/^ 2 (C)- Both the 
correlation length and the amplitude of G{r) decay to much smaller values when the line of 
MF instability is approached. The compressibility at low T* is very small for ( ^ 1/2 and 
very large otherwise. The small compressibility signals formation of regularly distributed 
clusters that repel each other upon compression of the system (when their separation be¬ 
comes shorter than the range of repulsion). The large compressibility signal the approach to 
the phase transitions between the periodic phase and the dilute gas or dense liquid phases 
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at T = 0. 

Our purpose is the verification if the mesoscopic theory in the self-consistent Gaussian 
approximation and with the further assumptions made in sec.4 is able to reproduce the 
above exact results. In order to calculate the correlation function and the /i*(C) and p*{C) 
lines, we have to solve Eqs. fl3^ - fl4T]l . To solve these equations we need the form of D{k), 
which depends on G{r) (see (155]) h In real space representation the correlation function 
G(o)(/c) = 1 /C 2 ^\k) for Cf\k) approximated by Eq. fl55D takes the form 

G^^\r) = cos(aor) -|- sin(aor) j, (46) 


where the correlation length is 


eo 



(47) 


and the wavenumber of the damped oscillations ao and the amplitude Aq are given by 


2 _ ^0 + ■\/fco + Cq/vq 

2 ’ 


(48) 


and 


1 


/l 2 = _ (' 40 ') 

AvoCoik^ + cq/vq) ■ 

Eq. (l46l) was obtained for continuum space. Our goal is to verify if our approximate theory 
reproduces the qualitative features of the exact results, therefore we assume the same form 
for the lattice model because of its simplicity. According to the ansatz (1551) . G{r) has the 
form given in Eq. fl45]) . but with the parameters Aq, replaced by the renormalized ones, 
Ar,^r,c^r respectively. The relation of the parameters Ar,^r,c^r with Cr,kr,Vr is analogous 
to the relation between Ao,^o,«o and co,ko,vo given above. The explicit form of D{k) can 
be easily obtained, and reads 


42 


^ “ ‘^Oir{k — 2ar) 


{k - 2ar)‘^ + 4^, 


-2 

r 


+ 


+ 2ar{k + 2ar) 4 ( q ;^ 4 -^^ ^) 


(fc + 2a,)2 + 4e 


-2 

r 


+ 


p + 4e-2 ■ 


(50) 


The above and Eqs. fl5U]) . fl45]) . fl55]) - fl57|) allow us to solve Eqs. fl55]) - flTT]) numerically. 

Before presenting the results we first verify if the assumption that C^\k) can be approx¬ 
imated by Eq. (l55|) is valid for this model. We Taylor expand V(k) given by Eq. fl44|) in terms 
of k^ about its minimum at k^ = k^. The form of V{k) lEq. (1541) 1 and the Taylor expansion 
truncated as in Eq. (155l) are shown in Fig.3 for J = 3 and J = 1/4. 
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FIG. 2: The lines of instability of the disordered phase in the Id lattice model of Ref.|l8l| in the MF 
approximation for the repulsion to attraction ratio J = 3 (solid line) and J = 1/4 (dashed line). 
Temperature is in units of the nearest-neighbor attraction (T* = ksT/ Ji). The volume fraction of 
particles C, is dimensionless. We study properties of the disordered phase inside the region where 
this phase is unstable in MF, i.e. below the shown lines. 



FIG. 3: V*{k) lEo. (11411 1 for J = 3 and J = 1/4 (solid lines) and the Taylor expansion leading to 
the approximate form of ()35p (dashed lines). The deeper minimum corresponds to J = 3. 
V*{k) is in units of the nearest-neighbor attraction Ji and k is in units of 1/cr, with a the particle 
diameter. 
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We can see that our approximation is good for k ^ kQ. However, the approximate formula 
overestimates and underestimates the effects of fluctuations with the wavelengths k < ko 
and k > ko respectively. Contrary to the assumption of a deep minimum (sec.4), in the case 
of J = 1/4 the minimum of V{k) is very shallow. 

In the second step we verify the ansatz (1381) . In our self-consistent Gaussian approxi¬ 
mation the inverse correlation function C{k) should satisfy both, Eq. (15H]) and Eg. (1^ . We 
compare Eqs. fl55]) and fIMD with each other, and with the ME version (15^ . In Fig.4a we 
show Eqs. (155]) . (15T]) and (I55|) for ^ = 0.24, T* = 0.1 and J = 3, and in Fig.4b for C = 0.41, 
T* = 0.1 and J = 1/4. Since —D{k) has minima aX k = 0,2ar (see (I50|) I and in the ansatz 
(l38|) there is a single minimum at /c = ko, the agreement between Eqs. (1381) and (I3T1) becomes 
worse for increasing IH 3 I, i.e. for increasing |^— 1/2|. Hence, the ansatz can be acceptable 
only for a limited range of </. 

For J = 1/4 it turns out that solutions of (I39p - (I41I) exist only in the central region 
of the ME instability of the disordered phase. For T* ^ 0.1 the solutions exist only for 
0.4 < C < 0 . 6 . As shown in Fig.4b, however, the agreement between (l38l) and (I3T1) is still 
not satisfactory for k < kr, except in the vicinity of ( = 1/2. Thus, the ansatz (l38|) is 
incorrect and either a different form of C{k) should be assumed, or numerical solution of 
dSH) is necessary. 

For J = 3 the agreement between the ansatz (1551) and the Eq. (15T]) is fair for a broader 
range of k. We shall consider 0.1 < C < 0.9, since our approximation (155|) becomes increas¬ 
ingly oversimplihed when small and large volume fractions are approached. We have verihed 
that the set of equations (159l) - (l4Tl) has solutions with > 0 below the line of ME instability 
shown in Fig.2. This indicates the lack of instability of the disordered phase, in agreement 
with exact results (see (155]) h 

Let us focus on the structure of the disordered phase. In Fig.5 we show the correlation 
length and the amplitude Ar of the correlation function as functions of ( for three different 
temperatures. The correlation function (see Eq. (1561) with renormalized parameters) is shown 
for C = 0.24 and two values of temperature, T* = 0.1 and T* = 0.7 in Fig. 6 . For low T* the 


correlation length is very large for ( 


/2, and rapidly decreases when If/ — 1/21 increases. 


in agreement with the results of Ref. 18j. On the quantitative level, however, the accuracy 
of the approximation decreases for decreasing T*. For T* ~ 0.5 we obtain semi-quantitative 
agreement between the approximate and exact results, but for T* ~ 0.1 the correlation length 
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FIG. 4: Comparison of the inverse correlation function C in MF (dash-dotted line), and as given by 
Eas. d^ and (I55|l . (dashed and solid line respectively). Top panel: J = 3, C = 0.24 and T* = 0.1. 
Bottom panel: J = 1/4, ^ = 0.41 and T* = 0.1. The parameters in (l38l) satisfy Eqs. p^ - (i4T]l . 
Note the negative values in MF, indicating the instability of the disordered phase, and the very 
small positive values in our theory. Note also that in the case of J = 1/4 the assumption (1381) is not 
satisfactory for k < kr- C is dimensionless and k is in units of l/u, with a the particle diameter. 


is significantly smaller and the amplitude is significantly larger than obtained in Ref. Q. 
Ar > 1 (Fig.6) is an artifact of our approximations. Note, however that in MF Aq ^ cxo for 
T* —)■ —V{kQ)/A 2 {C) from above, so the improvement in our theory is signihcant. 

In Fig.7 the wavenumber kr at the minimum of C{k) (maximum of the structure factor) 
is shown for J = 3 and T* = 0.4. It decreases slightly for increasing |(/ — 0.5[ , indicating 
increasing wavelength of inhomogeneities, in agreement with exact results Q . However, 


the magnitude of kr obtained in Ref. 


18| is slightly smaller. 


Let us focus on the effect of mesoscopic fluctuations on thermodynamic and mechanical 
properties in the Id model. The explicit expressions for the chemical potential and the EOS 
for fh given in Eq. (145 p have the form (see (13Op and 


p* = ■(/*(0)C + T* In 


C 


1-C 


+ T* 


and 


p* = -■(/*(0)C^ - T* ln(l -C)+T^ 


g(2C-l) 

2 ( 2 ( 1 - 0 ^ 

g(2C-l) 

2C(i-C)2' 


(51) 


(52) 
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FIG. 5: The correlation length (left panel) and the amplitude (right panel) as functions of 
the volume fraction ^ (dimensionless) for J = 3. From the top to the bottom line T* = 0.1,0.4,0.7. 
is in units of l/u, with a the particle diameter, and Aj. is dimensionless. 



FIG. 6: Correlation function for ( = 0.24 in the case of J = 3. T* = 0.1 (left) and T* = 0.7 (right). 
G(r) is dimensionless and r is in units of the particle diameter a. 


As discussed at the end of sec.4, we have neglected the last term in Eq. fl43H . 

When the fluctuation contribution (the last term in Eq. fl5T|) l is neglected, the slope of 
the line low T* is determined by E*(0). The system energy for ( = const, E(0), 

is positive for the strong repulsion to attraction ratio J = 3, and negative for J = 1/4 (see 
Fig.3), therefore in ME the slopes of the /i*(C) line are positive and negative for J = 3 and 
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FIG. 7: The wavenumber kr corresponding to the maximum of the structure factor (minimum of 
C{k)) as a function of the average volume fraction (dimensionless) for J = 3 and T* = 0.4. kr is 
in units of l/u, with a the particle diameter. 



FIG. 8: The reduced chemical potential as a function of C, (dimensionless) at T* = 0.1 for J = 3 
(left panel) and J = 1/4 (right panel). Solid and dashed lines show Eq. ()51li with and without 
the fluctuation contribution (last term) respectively. A narrow range of C, is shown for J = 1/4 
because beyond the shown interval there are no solutions of (I39p - ()4ip . is in units of Ji, with Ji 
denoting the nearest-neighbour attraction. 
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FIG. 9: The reduced chemical potential fj,*, Ea. (l30l) . (in units of Ji with Ji denoting the nearest- 
neighbor attraction) for J = 3 as a function of the volume fraction ( (dimensionless) deep inside 
the MF stability region of the periodic phase. Dashed, dash-dotted and solid lines correspond to 
T* = 0.7,0.4 and 0.1 respectively. 


J = 1/4 respectively. In the former case 3-particle clusters separated by 3 empty sites are 
energetically favourable and no phase separation occurs at T* = 0, whereas in the latter 
case a separation into dilute and dense phases occurs at T* = 0 for /r* = —0.75 |l8|. The 
MF instability for J = 1/4 (negative slope of the line) is associated with the MF 

phase separation that in exact results is absent for T* > 0. For T* > 0 exact results 181 


show positive slopes of /r*(C) in both cases, in agreements with the results of our theory 
(solid lines in Fig.8). From the exact results it follows that for J = 3 the slope at ( ~ 1/2 
is very large and increases for decreasing T*, whereas for J = 1/4 the slope at ( ~ 1/2 
is very small and decreases for decreasing T*. In Fig.8 we show the effects of mesoscopic 
fluctuations on fi*{C) for J = 3 and J=l/4atT* = 0.1. The very large difference between 
the slopes of the /i*(C) lines at (/ = 1/2 for J = 3 and J = 1/4 agrees with exact results. 
The compressibility at (/ = 1/2 is very small for J = 3 and very large for J = 1/4; the ratio 
between the compressibility is ~ 40. 

Let us describe the case of J = 3 in more detail. The slope of the lines for 

0.4 < C < 0.6 is almost independent of T*. Thus, our approximate theory does not reproduce 
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FIG. 10: Pressure p* (in units of the short-range attraction Ji per particle volume) as a function 
of the volume fraction (dimensionless) for J = 3. In the left panel Ea. (l52l) is shown for T* = 0.1 
(solid) T* = 0.4 (dash) and T* = 0.7 (dash-dotted line). In the right panel the effect of fluctuations 
is highlighted by comparison of the pressure with the MF result (dashed line) for T* = 0.1. 


the decreasing compressibility for decreasing temperature for ( ^ 1/2. On the other hand, 
when ( < 0.3 and ( > 0.7 the decreasing slopes of the lines with decreasing temperature 
are correctly predicted (see Fig.9). This behavior indicates that the compressibility increases 
to very large values at very small T* for small and large volume fractions. The increasing 
compressibility for decreasing T* signals the approach to the phase transitions that occur at 
T* = 0. The effect of the last term in Eq. flSU]) is clearly seen in Fig.8a. For J = 3 this term 
leads to a much smaller slope of the /i*(C) line for large- and small volume-fractions, and 
to a much larger slope for the volume fractions ~ 1/2, than in the absence of mesoscopic 
fluctuations. Note also that since /i* —)■ —cxo and /i* ^ cxo for /■ —)■ 0 and C ^ 1 respectively, 
three inflection points at the line /i*(C) must be present at low T*. This is a characteristic 
feature of the systems with strong mesoscopic inhomogeneities. 

Most of the above results agree qualitatively with the exact results in the Id lattice model. 
On the quantitative level, however, the accuracy of our predictions decreases for decreasing 
T*. The exact results for J = 3 show much smaller compressibility at ( = 1/2 for low T* 
than obtained in our approximation. Moreover, the very small compressibility at ( ~ 1/2 
increases to a very large value in a range of ( that is much more narrow than shown in Fig.8. 
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Finally, the isotherms /n*(C) intersect at 3 points: /i* = —2/3,2,14/3 |^, whereas in our 
approximation they are tangent to one another at fi* = 2. 

In Fig. 10 we present the EOS for J = 3. As in the case of the chemical potential, we 
obtain qualitative agreement with the exact results. In particular, at low T* the slope of 
the p*{C) line in the neighborhood of </ = 1/2 is much larger than for ( < 0.4 or ( > 0.6. 
When T* increases, the p*{C) line becomes smoother, in agreement with Ref.[^. However, 
the changes of the slope of the p*{C) line found in Ref. [l^ are much more pronounced at low 
T* than our theory predicts. Our isotherms do not intersect, in contrast to the isotherms 


obtained in Ref. 


18| . where they intersect in two points. 


VI. DISCUSSION 


We have combined DFT and field-theoretic methods in a coarse-grained theory for sys¬ 
tems with mesoscopic inhomogeneities [u, 3, 27]. Equations for the average volume 


fraction and the correlation function, fl28ll - fl2^ with fl25|l have been obtained in the self- 
consistent Gaussian approximation equivalent to the self-consistent 1-loop approximation in 
the field theory. Predictions of the theory for the disordered phase, fl5n]) - flM|) . have been 


verihed by the comparison with the exact results obtained in Ref. 


18| for a Id lattice model 


with hrst-neighbor attraction and third-neighbor repulsion. In order to simplify the calcu¬ 
lations, we have made further approximations concerning the inverse correlation function in 
MF, Eq. fl35l) . and beyond, Eq. fl38|) . and adopted for the correlation function in the lattice 
model the expression fl46]) obtained for a continuous system. Thus, we have in fact compared 
a continuous system with the lattice-gas form of the entropy, (I4S]) . with the exact results 
obtained for a lattice model. 

Despite all the simplifying assumptions, we have obtained a qualitative agreement for all 
the key features of the system with mesoscopic inhomogeneities, for the phase-space region 
where the ansatz fl38|) is a reasonable approximation. First of all, we predict stability of 
the disordered phase for T > 0. However, when the repulsion-to-attraction ratio is large, 
J > 1, the disordered phase is strongly inhomogeneous in this part of the phase-space region 
where MF predicts stability of the ordered periodic phase. In the ordered phase regularly 
distributed clusters are separated by voids. In the presence of mesoscopic fluctuations the 
long-range order is absent, but the inhomogeneity is reflected in the oscillatory decay of the 
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correlation function. The correlation length is very large and increases for decreasing 
T. Both the correlation length and the amplitude of the correlation function decay rapidly 
when the MF transition from the periodic to the homogeneous phase is approached (Fig.5 
for J = 3). All these features agree with exact results. 

The /i(C) and EOS isotherms also have the characteristic features observed in Ref. 
in the case of strong repulsion (J = 3). For low T there are 3 and 2 inflection points in 
the fi{() and p{() curves respectively. A very large compressibility for the volume fraction 
C ~ 0.3 becomes small for ( ^ 0.5, and very large again for ( ^ 0.7, in agreement with 
Ref.[^. The very small compressibility results from the repulsion between the clusters when 
the separation between them decreases upon compression. The very large compressibility 
accompany the approach to the phase transitions that occur at T = 0. The characteristic 
shape of the osmotic pressure can be used as an indication of inhomogeneities in experimental 
studies, when scattering experiments are not possible. We should note that in the case of 
weak repulsion, when the phase separation between dilute and dense phases occurs at T = 0, 
the compressibility at C ~ 0.5 is very large when T —)■ 0, in a striking contrast to the system 
with the strong repulsion between the particles, J = 3 (Fig.8). The only qualitative feature 
that is not correctly reproduced by our theory is the decrease of the compressibility at 
C ~ 1/2 upon a decrease of temperature for J = 3. 

On the quantitative level the accuracy of our results decreases for decreasing T. At 


very low T we obtain signihcantly smaller and larger amplitude than in Ref. 


18|. We 


should stress that the correlation function averaged over mesoscopic regions is not expected 
to be equal to the standard microscopic correlation function. Moreover, our numerous 
approximations applied to simplify the calculations influence the results on the quantitative 
level. Another source of discrepancy between the self-consistent Gaussian approximation 
and the exact results is the neglected effect of mesoscopic fluctuations on the higher-order 
correlation functions. Finally, the results could be systematically improved within the held- 
theoretic perturbation expansion beyond the self-consistent one-loop approximation. Still, 
the self-consistent Gaussian approximation gives results that agree qualitatively with the 
exact solutions, in contrast to MF, where a non existing phase transition and a divergent 
amplitude of the correlation function are obtained. 

We conclude that the disordering effect of mesoscopic fluctuations, neglected in MF, is 
overestimated in our theory, especially at low T. We expect signihcantly better results 
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without the simplifying assumptions, fl5^ . fl55]) and (SS]), but rather involved numerical 
computations are necessary to solve Eqs. fl5II]) - fl51|) . 

Our results show that most of the qualitative features of the disordered phase with strong 
inhomogeneities are correctly predicted even when very incomplete information about the 
interaction potential is taken into account (see the assumption fl35|) h Thus, the same proper¬ 
ties are obtained at this level of approximation for all interaction potentials which in Fourier 
representation have the same shape near the minimum. More generally, from our theory 
it follows that for C < Cc ( C > Cc)) where (c is the critical volume fraction, the chemical 
potential and pressure (see fl30|l and fl42ll i are smaller (larger) than predicted by MF. The 
fluctuation correction is the larger the stronger are the inhomogeneities, as measured by 
G(r,r). 

We hope that in future studies our self-consistent Gaussian approximation can be applied 
to systems with d > 1, where the exact results are not feasible. The present comparison 
with the exact results can help to interpret and critically analyze the results. 
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